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Abstract. We consider the real eigenfunctions of the Schrodinger operator on graphs, 
and count their nodal domains. The number of nodal domains fluctuates within an 
interval whose size equals the number of bonds B. For well connected graphs, with 
incommensurate bond lengths, the distribution of the number of nodal domains in the 
interval mentioned above approaches a Gaussian distribution in the limit when the 
number of vertices is large. The approach to this limit is not simple, and we discuss it 
in detail. At the same time we define a random wave model for graphs, and compare 
the predictions of this model with analytic and numerical computations. 

1. Introduction- the Schrodinger operator on graphs 

The structure of the nodal set of wave functions reflects the type of the underlying 
classical flow. This was suspected and discussed a long time ago, [TJ E3 El - and returned 
to the focus of current research once it was shown that not only the morphology, but 
the distribution of the number of nodal domains, is indicative of the nature of the 
underlying dynamics [S]. This was followed by several other studies of nodal statistics 
[HI 13 El EH EH- and their relation to the random waves ensemble [TT] . 
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Quantum graphs are excellent paradigms of quantum chaos |T2|, and in the present 
work we try to check to what extent the statistics of nodal domains in graphs follow 
the patterns observed in the study of wave functions of the Schrodinger operators in the 
typical systems (eg, billiards) where quantum chaos is often discussed. Graphs are one 
dimensional systems. Their complex features stem from two facts: i. their topology is 
different from a one-dimensional interval (except for starlike graphs they are generally 
not simply connected), ii. the corresponding "classical" dynamics is not deterministic. 
Because of the different topology, Sturm's oscillation theorem |T3] does not apply for 
graphs. We shall show however, that Courant's generalization of the oscillation theorem 
[Tl] to higher dimensions applies, but much more can be said about the problem. That 
is, the number of nodal domains of the n'th eigenfunction is generically bounded between 
n and n m i n and an explicit expression for n m i n will be given. 

The rest of this section will be devoted to the introduction of metric graphs and 
the corresponding Schrodinger operator. The nodal counting problem will be defined 
in section 121 and some general results will be presented. The distribution of the number 
of nodal domains will be discussed in section El This distribution will be calculated 
for star graphs in section 13.11 and the results of these computations will be used to 
derive the asymptotic distribution of the number of nodal domains in the limit of large 
graphs. Some results on counting domains on a bond in starlike graphs will be presented 
in section 13.1.21 Finally, in section 0] we shall introduce the random wave model for 
graphs. The mean and variance of the distribution of the number of nodal domains will 
be computed explicitly and compared with a few numerical results. 

A graph Q consist of V vertices connected by B bonds. The V x V connectivity 
matrix is defined by: 

Cij = number of bonds connecting the vertices i and j. (1) 

A graph is simple when for all i, j : Cy G [0, 1] (no parallel connections) and = 
(no loops). The valence of a vertex is = YlJ=i^i,j an d the number of bonds is 
B = | . =1 Ci t j. We denote the bond connecting the vertices i and j by b = 
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The notation and the letter b will be used whenever we refer to a bond without 
specifying a direction: b = [i,j] = To any vertex i we can assign the set of 

bonds which emanate from it: 

5« = {all bonds [i, k] : C i>k = 1}; = v t . (2) 

We assign the natural metric to the bonds. The position x of a point on the graph 
is determined by specifying on which bond b it is, and its distance x b from the vertex 
with the smaller index. The length of a bond is denoted by L b and, < x b < L b . 

The Schrodinger operator on Q consists of the one dimensional Laplacian on the 
bonds, which must be augmented by boundary conditions on the vertices to guarantee 
that the operator is self-adjoint. We derive the form of the boundary conditions here, 
since this way we can introduce several of the concepts and definitions to be used later on. 
Let x G Q and *&(x) a real valued and continuous function on Q, so that ty(x) = i^ b (x b ) 
for x G b, and < x b < L b . The functions ip b (x b ) are real valued, bounded with piecewise 
continuous first derivatives. The set of functions which fulfill these conditions will 
be denoted by V and they are the domain of the (positive definite) quadratic form 

Q[¥] = jf dx(Vn*)) 2 = E dx b (^j 2 . (3) 
The unique self-adjoint extension for the Schrodinger operator, H, is determined by the 
Euler - Lagrange extremum principle. The domain of H, T>h consists of functions in V, 
with twice differentiable if) b (x b ), which satisfy the boundary conditions 

0, (4) 



\/i = l,...,V : £ 



n b (i 



dx b 

beSM 

where the derivatives are computed at the common vertex, and n b (i) takes the value 1 
or —1 if the vertex i is approached by taking x b to or L b , respectively. These boundary 
conditions are referred to as the Neumann boundary conditions. In the following we 
shall denote by <pi the value of \& at the vertex i. 

The spectrum of the Schrodinger operator H is discrete, non negative and 
unbounded. It is computed by solving 

-^^ = *2* W <*), V* (5) 
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subject to the boundary conditions (J1J). The resulting eigenvalues are denoted by by k^, 
and they are ordered so that k n < k m if n < m. 

For later use we quote the following property. Let V n denote the subspace of 
functions in T> which are orthogonal to the first n — 1 eigenfunctions of H . Then, for 
any non zero $ G T> n 

Q[$] > k\ f dx($ 2 (x)) . (6) 
Jg 

Equality holds if and only if $ is the n'th eigenfunction of H . 

It is convenient to present the solutions of © on the bond b = {i < j) as 

i>b( x b) = . -j T {4>j sin kx b + 4>i sin k(L b - x b )) . (7) 
sin kL b 

The spectrum is computed by substituting ((7j) in (jlj), which results in a set of linear and 
homogeneous equations for the fa with k dependent coefficients hij(k). The spectrum 
is obtained as the solutions of the equation £(k) = det h(k) = 0. As will be explained 
shortly, we shall assume that the lengths L b are rationally independent, so that ((k) 
is an almost periodic function of k. We shall also restrict our attention to simple and 
connected graphs, and to avoid lengthy discussions of special cases, will assume that 
the valences Vi > 3 at all the vertices (exceptions will be stated explicitly). 

2. Nodal domains on graphs 

Nodal domains are connected components of Q where the wave-function has a constant 
sign. One cannot exclude the possibility that eigenfunctions of the Schrodinger operator 
vanish identically on one or several bonds. This is often the case if the bond lengths are 
rationally dependent. As an example, consider three vertices which are connected by 
bonds which form a triangle. If the lengths of each of the bonds are integer multiples of 
L t , then there exist eigenfunctions with eigenvalues k t , n = n|^,for any integer n which 
vanish on all the other bonds ^j- To exclude such cases, we shall discuss graphs with 
lengths which are rationally independent (incommensurate). 
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Rational independence is not sufficient to remove completely the possibility that 
wave functions vanish along one or several bonds. To construct an example, take any 
graph and choose a wave function which has a few nodal points on it. Connect the nodal 
points by bonds and take their length such that the new graph has incommensurate 
lengths. The Schrodinger operator for the newly constructed graph has the same 
eigenvalue, and a wave function which vanishes identically on all the added bonds. 
This construction is quite general, but at most, it can bring about a negligible number 
of such wave functions. The reason for this is as follows. The n'th wave function on the 
bond vanishes if and only if both (jJf 1 ^ and vanish. The vectors (f)^ are the 
null vector of the quasi-periodic matrix hij(k n ), and as k n goes over the spectrum, they 
cover the sphere ergodically. Thus, the probability that several components are exactly 
is vanishingly small. In the sequel we shall ignore these non-generic cases, but bear in 
mind, however, that their presence cannot be completely excluded. 
The nodal domains on graphs are divided into two types: 

i. interior domain - A domain which is restricted to a single bond , and whose length 
is exactly |. 

ii. vertex domain - A domain which includes a vertex, and extends to the bonds which 
emanate from it. 

There are V vertex domains, and their length Aj can take any value in the range 
ViT > Aj > 0. Denoting the length of the graph by £ = Ylb=i -^ fe > we obtain the following 
expression for the number of nodal domains, 



7T 



= 1 



Note that the second term above is an integer, and that this expression is correct for the 
generic case where the wave-function does not vanish along entire bonds. v n is bounded 
in the interval 

b£ + V > v n > ^ + V - 2B . (9) 

7T 7T 

In the limit n — > oo, — > I. Hence, in this limit, — — > 1. This observation 
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stands intermediately between Sturm's oscillation theorem = 1), and Pleijel's result 
that lim^ is strictly smaller than 1 for the eigenfunctions of the Dirichlet Laplacian for 
domains in M 2 |16j . 

An alternative expression for the number of nodal domains provides a sharper 
bound on the range of variation of v n . Denoting the number of nodal points on the 
bond b = by jin\ we have 



u (b) 

H'n 



- (l - (-l)^lsign[^]sign[0,]) (10) 

where [x] stands for the largest integer which is smaller than x, and fa, <pj are the values 
of the eigenfunction at the vertices i,j respectively. Thus, 

B 

Vn = J2»n- B + V - ( U ) 

6=1 

The allowed range for v n is now 

E(— i + v * v « > D— ] + v - B ■ ^ 

6=1 6=1 



The estimates from above can be sharpened by Courant's law [M] adapted for the 
present problem, which we shall now state and prove. 

Theorem: Let Q be a simple, connected graph. Let k\ be the n'th eigenvalue of 
the Schrddinger operator H defined above and let ^ n (x) be the corresponding real 
eigenfunction. Then, the number of nodal domains v n of^ n [x) is bounded from above 
by n, and this bound is optimal. 

The proof follows the method used in [16] . Assume that v n > n. Denote by 7/, 
the nodal domains on Q, so that U£i 7' = @- Construct n functions Ui(x) 6 D in the 
following way: 

{fy n (x) if x G 7 m 

(13) 
otherwise 

It is always possible to find n real constants a m such that U(x) = Ylm=i a mU m (x) 
is orthogonal to the first n — 1 eigenfunctions of H. Hence U(x) G T> n . A simple 
computation shows that 

Q[U] =%[ U\x) . (14) 
JQ 
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However, U(x) is not an eigenfunction, hence the above equality is in contradiction with 
the strong inequality imposed by (|Bj). Thus, the assumption that v n > n is false. □ 

In the next sections we shall try to determine how the v n are distributed within 
their allowed range. We shall start by solving a simpler problem, which pertains to 
the family of star graphs. A similar derivation for more complicated graphs is beyond 
our present ability. However, assuming that in the limit of large graphs, the lengths 
of vertex nodal domains are independent, we shall be able to deduce an approximate 
expression for the distribution of v n in this limit. 

3. Nodal domain statistics on graphs 

In the previous section we observed that ^ — > 1 as n — > oo. Hence, there is no point to 
use the statistics proposed in [5] for graphs. Rather, we shall discuss the distribution of 
the quantity 

8v n = v n -n (15) 

which can vary in the interval + V — IB — n< b~v n < 0. Let X n = ^ YlJ=i ^-t denote 
the sum of the lengths of the vertex domains measured in units of half the wave length. 
Following (jHJ) we find 

8v n 



k n C 1 

H n 



+ V- l --\ n = SN(k n ) + V - l - - \ n . (16) 



7T 2 

The expression in the square brackets above is the deviation 5N(k n ) of the mean spectral 
counting function from its actual value. Thus, the fluctuations in the number 
of nodal domains stem from two sources: the spectral counting fluctuations, and the 
fluctuations in the lengths of the vertex domains X n , whose distribution we shall denote 
by 

P(X) = (S(X-X n )) , (17) 

where (■ • •) indicates average over a spectral interval of size Ah, with eigenvalues 
on average. In general, and especially for graphs with small B, the two contributions 
are probably correlated. For large graphs, however, such correlations are expected to be 
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much weaker. This were the case if the spectral and the eigenvector distributions are 
independent, like in the relevant random matrix ensemble (GOE). We are not able 
to prove this statement, and we assume that in the limit of large graphs the two 
contributions can be treated independently. The quantities of interest here are 

M = ((^ + \- n ))+ V -\-(*») ( 18 ) 

and 

(A^ 2 ) - U {^f + \ - nj \ + (AA 2 ) . (19) 

The contribution of the spectral fluctuations to the mean (Su), vanishes O (/□:)• A n is 
bounded to < A n < IB. In the sequel, we shall provide evidence in support of the 
natural expectation that the A distribution is symmetric around the point B, thus 

(A) = B . (20) 

Hence, 

(v n -n) = -(B-V + ±). (21) 

This result is consistent with v n — n < since we assumed Vi > 2 at all the vertices. 

Turning to the variances, the contribution from the spectral counting function for 
general systems, and for graphs in particular, was studied previously by various authors. 
We show in Appendix (A) that 

<W(W> =f fl + of^Y) , (22) 



6 V V B 

The main term in (j22j) is a universal bound which is valid for all incommensurate 
Neumann graphs. The error estimate is valid for well connected graphs, where the 
spectral statistics is known to follow the random matrix predictions. The rest of this 
and the following sections will deal with the distribution of the total size of the vertex 
domains A n . 

The distribution P(A) for a finite small graph shows distinctive features as can be 
seen in figure (0) where we plotted the numerically obtained distribution for a fully 
connected graph with V = 4 vertices and B = 6 bonds (the tetrahedron) and compare 
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P(AX) 



P(AA) 




15 -10 -5 

AA = A - (A) 



Figure 1. Distribution of the deviation A A of the total length of vertex domains 
A from its mean value (A) = B for fully connected graphs with V — 4,5,8 and 10 
vertices (in units of half the wave length). The first 10 7 eigenfunctions have been used 
for the numerically obtained full line. The thin line is a Gaussian of variance B/6 



it to larger graphs. A bell shaped function is obtained for P(\) of larger graphs which 
is quite well approximated by a Gaussian of variance ^. 

Since A n is bounded in the interval (0, 2B) its variance cannot grow faster than B 

(AA 2 ) = /3(G)B , (23) 

where (3(G) < 1. We shall compute (3(G) below for particular models, and show that in 
general, for large graphs (3(G) = \- 

A quantity which might be of some interest in the present context is the length 
(again in units of half the wave length) Xn'^ of the intersection of a vertex domain at 
a given vertex % with the single bond b = Generally, Xn'^ ^ Xn'^ but they are 
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related by 



7T 



.k n L 



Mil _u v(*,i) -i_ vO'.*) 



\ (-i)i^i 

sign^] sign^- 



(24) 



2 2 
The total length of all vertex domains is 

a, »,,(C^ • \;/-) . (25) 

The distribution 

^ W) (X) = <*(X-XS J) )> (26) 

is thus connected to nodal counting on a single bond. Due to (strong) correlations 
between the Xn for different % and j ()26|) is less useful than (JT7j) or nodal counting on 
a complete graph. 

5.1. Nodal domain statistics on star graphs 

In a star graph all bonds emanate starlike from one central vertex % = 0. Each bond 
hi — [0, i] (i = 1, . . . , S) connects the central vertex to one peripheral vertex i (see figure 
12} . The bond 6; has the length Lj and all lengths are chosen incommensurate. The 
variable Xj measures the distance from the center on bond bi such that < Xi < Lj and 
Xj = Lj at the peripheral vertex z. Though the number of vertices is V = B + 1 only 
the central vertex fulfills vq > 3 (if i? > 3). The peripheral vertices have valence v$ = 1 
and instead of Neumann boundary conditions we will use Dirichlet boundary conditions 
<f>i = (i = 1, . . . , B) there. The wave function on the bond bi follows from f7j) 

i>i{xi)= . ° sink {Li - (27) 

Sill rZJ-Ji 

where <po is the value of the wave function on the central vertex. Current conservation 
at the center leads to the quantization condition 

B 

f B {k n )^J2 cot k ^ = ° ( 28 ) 
for the n'th eigenvalue k n of the star graph. 
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11 




Figure 2. A star graph with B = 6 bonds emanating from the central vertex 0. 

Since the peripheral vertices are nodal points, equations (JHJ) and (fTHjl for the number 
of nodal domains have to be modified to 

^ = l + E[^] = l + -£-A n (29) 

i=l 

where A n = — Xl^LiI^r 21 ] * s ~ times the length of the nodal domain that contains 
the central vertex. Obviously < X n < B and v n is bounded by 

l + h£> Vn >i- B + ^. (30) 

7T 7T 

3.1.1. The central nodal domains 

As discussed above, nodal counting is partly determined by spectral fluctuations and 
partly by the distribution ()17|) of the length A of the central nodal domain. We shall 
consider here only the distribution of the lengths of the central vertex domain, 

P{\) ={8{X-K)) 



TV 

lim / dk 







dk 

where \{k) = J2f=i - I^ 1 !)- Fr o m (El we have 



a/. , , (31) 
5(/ s (fc))5(A-A(A;)) 



d/s (A;) 



dA; ' sin 2 kL 

i=i 



Let 



X, = ^ - 1^1 (33) 

7T 7T 

be the (rescaled) length of the intersection of the central nodal domain with the i'th 
bond (X(k) = Y^i=iXi)- Obviously, < Xi < 1 an d since the length Li are assumed 
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incommensurate, k creates an ergodic flow on the £>-torus spanned by the x% EI]- Thus, 
the spectral integral over k in (pTTJ) may be replaced by an integral over the S-torus 
variables Xi- This leads to 



X* • (34) 



Replacing the two ^-functions by their Fourier representation, the distribution takes the 
form 

i roc roo 

P(X) = — dr, d£ G(r), O^Giv, (35) 

J — oo J — oo 



where 



G(V->0 = ~ I da cos I 77 tan a + — a ) (36) 

7T In \ 7T 



and 



6(77,0 = (l-J^Gfo,0 

= 2 cos ^ 5(77) - - G(t7, P- . 

2 7T 



(37) 



The last line shows, that G(r),£) is a distribution where denotes the Cauchy's 
principal value. The integral (}3l)j) can be solved explicitly (see [221, (3.718)) in terms of 
Whittaker functions W^ u (x) 

W_±_ ! (277) W X ! (-27/) 

Gfafl = 8(17) "r 3 ^ + e(-Ty) 7' g ^ • (38) 

Here G(x) is Heaviside's step function. The appearance of Whittaker functions can 
also be seen from (|37j) - for 77 7^ the right hand sides reduce to (1 — A) G(r],£) = 
— ^ G(rj,£), a special case of Whittaker's differential equation [221 • Using the last line of 
equation (|3*7j) the distribution can be written as a sum of two terms P(X) = Ps(X)+Pp(X) 
where 

P S (X) = — d£/ d7 / G(r / ,0 B - 1 ^)cosicos(e(A--)) 

27T J.oo J.oo 2 2 



7T 



00 / • \ B—l 



2 /"°° , /sinz . 

dz I cosz cos(z(2A — B)) 



(39) 

— E ^ Z ^(0(f-HA-f|)(f-HA-||) s - 2 + 

0<«<^ 



0(f-l-^|A-4|)(f-l^|A-4|) B - 2 ) 
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and 



1 f°° R f°° 1 

= -4^1 dacos(e(A--))y d^p-c^o 



2 

^2 



7T 



/■oo 

/ dz2cos(^(2A - B))x 
Jo 




(40) 



TV ' 2 



r(i - * 

P(A) is symmetric in A — Hence (A) = -f . 

For large star graphs (£> ^> 1), -P(A) is dominated by Ps(X) (see figure This 
observation is supported by the fact that Pp(A) does not contribute to the integrated 
probability distribution, J^dAP^A) = while J^dX Ps(X) = 1. The dominance 
of Ps(X) can be further supported by computing the variance of the exact distribution 
P(A) and of P^(A). The exact variance is evaluated by going back to (|34jl : 

(AA 2 ) =|dA (A-f) 2 p(A) 

J sin vrxi V tt J V 2 tT 



(41) 



P + 2 4 p tanz7r 
dz 2; arctan ■ 



12 7T J P - 1 

^ + 2 + O^B' 1 -") , (p>0). 



12 

Using (}3T?j) we reproduce the leading terms in the exact variance: 

/ dA ( A _g\ AH ^. (42) 

The fact that Ps{\) approaches P(A) for large star graphs is very important in the 
present context. First, it provides an analytic expression, which for large B tends to a 
Gaussian with a variance given by ()42j) . Second, when we consider general large graphs, 
the size of the vertex domains become statistically independent, and their distribution 
can be approximated by a Gaussian whose variance is 

((AA) 2 ) «-^5> = f , hence (3(G) = I (43) 

where here B stands for the number of bonds on the general graph. This result is 
consistent with the numerical data shown in figure ^ 
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Combining the two estimates for the variances of the spectral fluctuations (J22j) and 
the nodal domain fluctuations (j43jl . we obtain the leading term for the variance of the 
number of nodal domains: 

((A^) 2 ) = f . (44) 

This estimate holds in the limit of large graphs. In the next section we shall show that 
the random wave model for the graph provides the same answer for the variance of the 
nodal domain distribution. 

3.1.2. Nodal domains on a single bond 

In a star graph the number of nodal domains on the bond bi — [0, i] is 

I ] = Xi ( 45 ) 

7T 7T 

where Xi is the length of the intersection of the central nodal domain with the bond bi 
as in equation (J3*3*j) above. Note, that there are no vertex domains on the peripheral 
vertices and equation (|24j) has to be modified. Following the preceding section we define 

ii if. 

S(f B (k))6( X - X.W) 

u (46) 



7T f 

P {i \x) = lim -tt-t / dA; 



dfB (k) 



dk 



71 



C sin 2 irxj 

With similar techniques as in the previous section this integral can be solved explicitly 

P ( Hx) = ^ + ± j= U2 B ~\ 2 • (47) 
L L (B — ly sm x n + cos X 71 

For large star graphs B ^> 1 where each bond length is of similar size, one has ^ ~ -g, 
and the distribution becomes uniform (with two singular points at x — and x = 1). If 
one bond length Lj exceeds the other bond lengths such that L.; ^> C—Li the distribution 
becomes 

pW W w m W* B *~\ ^ (48) 

[B — l) z sm x n + cos X n 
which for large B is peaked at x = and x = n - 
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Figure 3. Numerically obtained distributions of the length of the central nodal domain 
in a star graph with B = 3, 4, 5, 50, 100 and 200 bonds (histograms - full lines) - 
the first 2 • 10 6 eigenfunctions have been used for each graph. The dashed line gives 
P^(A) = P(A) — -Pp(A) and the thin line is a Gaussian with variance y| - for large B 
the Gaussian is indistinguishable from the -P(s(A) and the numerical distribution slowly 
converges to the Gaussian. 



4. Random waves on graphs 

Since Berry's seminal work ^Tj random waves have been a paradigm for chaotic wave 
functions. Recently they have been used extensively in the investigation of the nodal 
structure in chaotic wave functions 13 El El HO] ■ 

In this section we will introduce random waves on graphs. Any ensemble of random 
waves should solve Schrodinger's equation on the bonds and be continuous at the 
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vertices. Thus, any set of values {(f) j} (j = 1,...,V) for the wave function on the 
vertices determines a wave function on the graph which solves (JJJ). However, these waves 
do not fulfill the correct boundary conditions (current conservation) on the graph. The 
ensemble of random waves on a graph is therefore defined in terms of the probability 
distribution of the <pj. 

Before we define the appropriate ensemble for more general graphs it will be 
instructive to discuss star graphs. If we want to compare random waves with the star 
graph results of the previous chapter we have to keep the Dirichlet boundary conditions 
at the peripheral vertices and only relax the the current conservation condition at the 
center. Then tfji(xi) = si ^ L . sin fc(L, — Xj) is the random wave on the z'th bond and (fi is 
the only random parameter for fixed k. Obviously the position of nodal points does not 
depend on <p . Thus, we do not need the distribution P(<fio) for the discussion of nodal 
domain statistics. Let us now rewrite equation (|29j) for the number of nodal domains 

u(k) =i + y d Xt v 8(xi-j) 



i=l j=-oo 



kC B -J^. sin 2tt jLik 
= 1 + t^- ■ 

i=l j=l J 

Averaging over a fc-interval reveals that the mean is [u) — 1 + — — ^ ~~ consistent with 
the result for the eigenfunctions in the previous chapter. For the variance we get 

cos2irk(jLi—j'L i /)—cos2iTk(jLi-\-j'L i i) 



(^) =1 EE 



kn "T 1 2Akjj'n 2 



which again coincides with the result for the eigenfunctions of star graphs for large 
B. Note, that for incommensurate bond lengths there are no correlations between 
the contributions from single bonds. Thus the number of nodal domains is a sum 
of independent quantities of finite variance each. The central limit theorem leads to 
Gaussian statistics for large B. 

To define an appropriate ensemble for more general graphs we will be guided by 
wave functions that do fulfill current conservation (we will again assume that each vertex 
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has a valence Vi > 3). The quantization condition for a graph has the form det hij(k) = 
where hij(k) is a real symmetric matrix of dimension V x V ^2]- If the quantization 
condition is fulfilled for k, the eigenvector for the zero eigenvalue is a set of vertex values 
{4>j} that determines the eigenfunction. For incommensurate lengths hij(k) is a quasi- 
periodic function of k such that the matrices hij(k) are expected to be typical members 
of the Gaussian orthogonal ensemble (GOE) in random-matrix theory. Since for that 
ensemble eigenvectors have uncorrelated components, we will assume for the ensemble 
of random waves on the graph that the cj)j are independent Gaussian variables of equal 
variance. From equation (jlOjl and (jllj) we see that the number of nodal domains only 
depends on the signs of 4>j- Let o~j = sign[</> 3 -], then Oj = ±1 with equal probability. The 
number of nodal domains can now be rewritten as 



kC _ v — -\ „ / v — > sin2irTnkLa 1 , . r kL ij i \ , . 

K* = — + V - B + E C m E - ^(-l) 1 — 1 ^ -(51) 

«j':*<i \m=l / 

Averaging over a ^-interval of length Ak centered at k and over o~j gives the mean 
number of nodal domains 

(u) = — + V - B + OiAk' 1 ) . (52) 

7T 

The variance of the number of nodal domains is purely due to the sum over i,j in ([51)1 . 
To leading order, the sum over sines and the sum over the signs Oi give independent 
contributions to the variance. We have already calculated the first part (^ 2 ) S m = y| 
above in our discussion of random waves on star graphs. The fluctuations due to the 
signs are stronger, and they provide to the variance a term (z/ 2 ) CT = j. Hence the 
random wave model predicts the variance 

<^ 2 > = f + OiAk- 1 ) . (53) 



This result reproduces the estimate (|44|). which was derived under very different 
assumptions. Finally we would like to note that bonds do not contribute independently 
to the number of nodal domains. However, these correlations are not expressed in the 
variance. They do contribute to the higher moments. 
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6. Appendix A 

In this appendix we justify the bound (}22|) on the variance of the spectral counting 
function for graphs. 

The starting point is the expression of the spectral counting function N(k) as a 
sum of its mean value (Weyl's law) and an oscillatory part, 

kC 1 

N (k) = — + - + SN(k), (54) 

7T 2 



and the oscillatory part is given by [12] 

OTW = I Im f;!l(M0r. (55) 

m=l 

iSb(^w) is the bond scattering matrix defined in [T2*] . tr (iSb(^))"\ is a sum of 
contributions from all the m— periodic orbits on the graph. Each contribution is endowed 
with a phase proportional to kl™, where 1™ is the length of the orbit, and p is the 
summation index. Because of the rational independence of the bond lengths, 

(tr {S B (k)) m tx (S B (k)) n ) = 2BK m S m>n + <D(Ak~ l ) (56) 

Where K m is the spectral form factor associated with the graph. The mere fact that the 
bond lengths of the graph are not commensurate, is enough to guarantee that K m < 1. 

z7r L — ' m z o 

m=l 

A sharper estimate can be given for well connected graphs, where numerical and analytic 
results O EI] show that the K m follow the predictions of random matrix theory for 
the circular orthogonal ensemble (COE). We can use the known functional dependence 
of K m on m and B and show that 

((^)) 2 ) = f(l + ^^)- (58) 
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